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Abstract 

I In this paper, we study the rate of convergence in periodic homogenization of scalar ordinary differential equations. We 
' provide a quantitative error estimate between the solutions of a first-order ordinary differential equation with rapidly oscil- 
OO lating coefficients and the limiting homogenized solution. As an application of our result, we obtain an error estimate for 
the solution of some particular linear transport equations. 

<\ 

^ 1 Introduction 

a 

g 1.1 Homogenization of an ODE 

In this paper, we consider the solutions of the following first-order ordinary differential equation: 



> 



X 




t > 0, 

(i.i) 

[ u'(0) = -Uo, 

where e > 0, stands for or equivalently dtu'^, and uq is a real number. We are interested in the rate of 



convergence of the solution u'^ to its limit in the framework of periodic homogenization. We employ the following 



m 

' assumptions on the function / 

o 

> 



• (Al) Regularity: the function / : ^ K is a bounded Lipschitz continuous function with a — Lip{f) its 
Lipschitz constant, and f3 = ||/||loc(r4); 

• (A2) Periodicity: for any {v,T,u,t) G M^, we have: 

f{v + l,T + k,u,t) — f{v,T,u,t) for any {I , k) € 7? ; 

• (A3) Monotonicity: for any (w,r, i) G R^, 

the function u i — > f(v, r, u, t) is non-increasing. 



Let us make short comments on these assumptions. Remark that assumption (Al) ensures the existence and 
uniqueness of the solution u'^ of (|l.ip via the Cauchy-Lipschitz theorem. Moreover, the assumed boundedness of / 
is not a restrictive condition while we work on any finite time interval [0,T]. The monotonicity assumption (A3) 
may seem unnecessary at a first glance, but will be indeed useful to guarantee the uniqueness of the solution to the 
homogenized equation (see Proposition II. 4p . Moreover, assumption (A3) will play a crucial role to establish the 
rate of convergence of u'^ to its limit u° (see for instance Section S]). 

In order to define the homogenized equation, we will use the following proposition: 
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Proposition 1.1 (Definition and properties of the effective slope 0f). Fix {u,t) G M^. Then there exists 
I G K. such that for any initial data Uq G M. the solution v G C"'^([0, oo); M) of the following ordinary differential 
equation: 

j Vr ^ f{v,T,u,t), r > 0, 
I V{0) = Mo, 

satisfies 

as T^oo. (1.2) 

T 

Let us set the effective slope: 

IM^I. (1.3) 

Then the following holds: 



I 0f : K'^ — !■ R is continuous. 

{ (1-4) 

I For any t>0, the map u ^ 0f{u,t) is non-increasing. 

Let us mention that, for some specific functions /, explicit formulas for 0/ can be obtained (see for instance [23] . 
and the examples below). 

Example 1.2 For f{v, t, u, t) = —u + cos(27ru), we have 0f{u, t) = -M+cos(27r'») ) f'^^ ^ ^ ^' "^""^ 

0/(m, t) ~ —c\Ju — 1 as ti^l^, 

/or some constant c > 0. 

Example 1.3 for /(u, r, u, t) = — u + | sin(27ru)|, we have, for some constant c > 0: 

0/(u,O~TT— m "■^ w->0-. (1.5) 

Example 11.21 shows in particular that even for analytic /, the function 0/ could be non-Lipschitz. Example 11.31 
shows a case where 0/ is not Holder continuous when / is Lipschitz continuous. The proof of Example 11.31 will 
be given in the Appendix. At this stage, we can write the homogenized equation associated to equation (jl.ip as 
follows: _ 

ru? = /(w°,t), t>o, 

1 u"{0) ^uo. 



Even if 0/ may not be Lipschitz continuous in u, we can show the existence and uniqueness of the solution of p.6p . 
taking advantage of the monotonicity of 0/(u°,t) in Indeed, we have: 

Proposition 1.4 (Existence and uniqueness). Under assumption on 0f , there exists a unique solution 

u° G Ci([0,oo);M) of ra . 

It is worth noticing that assumption p.4p satisfied by 0/ in the homogenized equation is our motivation to make 
assumption (A3) on /. We can now state our main result: 

Theorem 1.5 (Error estimate for ODEs). Under assumptions (Al)-(A2)-(A3), if u'^ is the solution of 

and is the solution of the homogenized equation \1.6\) . then for every C > 0, e > 0, and every T > Ce\ loge|, we 

have the following estimate: 

K--°IU-(o,T)<^^, (1.7) 
where c > is a positive constant only depending on C and on a, (3 defined in assumption (Al). 
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Such a result for a general monotone system of ODEs seems to be completely open. The above estimate in ji^^^ 
is in fact related to the behavior of 0/ in Example 11.31 which is the worst possible regularity of 0/. Moreover, it is 
possible to show that under the condition T > Ce\ loge|, inequality (|1.7[) is sharp, see the following example whose 
proof will be given in the Appendix: 

Example 1.6 Let f{v, r, u, t) = g(v + t) — 1 with a 1-periodic function g satisfying 

g{w)^\w~l/2\ for we [0,1]. (1.8) 

In this case 0f{u,t) = —1. Let us choose the initial data uq = 0. Then for any 6 > 0, we have the following 
estimate between the solution u'^ to and to \L6\l : 

u\t)-u^{t)^,^^^^ for t^5e\\oge\. 

Remark 1.7 It is worth mentioning that assumption (Al) could be replaced by the weaker assumption: 

• (Al)' Regularity: the function / : R is a bounded continuous function such that for every r G R, the 

function f{.,T, ., .) is Lipschitz continuous. 



1.2 Application to the homogenization of linear transport equations 

For a; = (xi, X2) e R^, let us consider a vector field a'^ = (af , a|) defined as follows 

ai{xi,X2) = -f y—,—,xi,X2j ^^^^ 

4(2:1, 0:2) = 1, 

with a function / satisfying (A1)-(A2)-(A3). We consider the viscosity solution V'^{t,x) of the following linear 
transport equation: 

f + a' ■ ^0 on (0, 00) x R^ 
[V'{0,x) = Vo{x) on R^ 

where Vo : R^ — + M is a Lipschitz continuous function. The existence and uniqueness of a viscosity solution 
of (jl.lOp is ensured since a*^ g PF^'°°(M^) and Vq is Lipschitz continuous (see for instance [3]). The expected 
homogenized equation associated to (|1.10p is: 

f VP + a- VF° = on (0, 00) x R^ 

<^ (1.11) 

\v°{0,x) ^Vn{x) on R^, 

with the vector field 0a = (001,002) defined as: 

J 0ai(xi,a;2) = -0/(a;i,X2), with 0/ given by (fO)l . 
I 002(^1, 2:2) = 1. 

As a consequence of Theorem II .51 we will show in Section [B] the following result: 



(1.12) 



Theorem 1.8 (Error estimate for linear transport equations). Under the previous assumptions, there exists 
a Lipschitz continuous function which is a viscosity solution of such that for any C > 0, e > 0, and 

T > Ce\ loge|, the solution of U.10\) satisfies: 

||y^-y°|U»(K.,(o,T))<^ (1.13) 



where c' = cLip{Va), and c > is the constant given in Theorem \1.5\ 
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Choosing the initial condition Vq{x) = xi, we can easily deduce from Example 11.61 that inequality p. 131) is also 
sharp for T > Ce\ loge|. Here, in this application, the vector field a*^ is quite special. The interested reader could be 
referred to [14] for some other examples of vector fields in 2D, where a homogenization result is presented without 
any rate of convergence. In [27j , the author gives some non-explicit error estimates for linear transport equations in 
the particular case of periodic vector field a*^. However, these error estimates obtained in [27] may depend strongly 
on the irrationality of the rotation number loq associated to the vector field a'^ (where wq is nothing else than —0/ 
in our application). On the contrary, estimate (|1.13p only depends on some bounds of the data of the problem, and 
are completely uniform with respect to the rotation number. 

Remark that when /(w, t, u, t) is independent of u and t, we have a much better estimate: 

Theorem 1.9 (Better error estimate). Under the assumptions of Theorem \ 1.81 if f{v,T,u,t) is independent 
of u and t, then we have: 

I1^^'-i^"IIl-(r^x(o,t)) <c"e, vr>o, 

where c" ~ ^Lip(Vo), with ^ (given in Provosition \2. 1\) only depends on [3 defined in assumption (Al). 
The proof of Theorem 1 1.91 will also be given in Section [S] 

1.3 Brief review of the literature 

The pioneering work (via the theory of viscosity solutions) to periodic homogenization was established in [17] . 
Starting from [17j . the homogenization theory for Hamilton- Jacobi equations has received a considerable interest. 
There is a huge literature that we cannot cite in details, but the interested reader can for instance see [TJ [SI HI [TTl 
[T51 [TSl [TB] and the references therein. Another aspect concerning homogenization of SDEs (stochastic differential 
equations) has also been studied by several authors (see for instance [^IT51 [51[^V These problems are related to 
our problem when the SDK reduces to an ODE. 

To our knowledge, the question of estimating the rate of convergence in homogenization of PDEs has not been 
widely tackled up elsewhere in the literature. We can cite [7^ for several error estimates concerning the rate of 
convergence of the approximation scheme to the effective Hamiltonian. We can also cite the work in [8| about the 
rate of convergence in periodic homogenization of first-order stationary Hamilton- Jacobi equations, where an error 
estimate in e^^^ is obtained for Hamilton-Jacobi equations with Lipschitz effective Hamiltonian. 

For the problems of homogenization of ODEs, we refer the reader to [23l [24] . We also refer the reader to 
[2 [9l [TOl [T9[ [211 [221 [26] for problems on homogenization of nonlinear first-order ODEs and/or the associated linear 
transport equations. As mentioned above, we refer the reader to [27| for some other error estimates for linear 
transport equations. 

1.4 Organization of the paper 

The paper is organized as follows. In Section [21 we present the proof of an crgodicity result fProposition l2.ip that 
defines 0f = I. We also present the proofs of Propositions ll.il and ll.4l In Section [S] we give a result of stability of 

1 under additive perturbation (Proposition 13 . ip . A basic error estimate (Proposition [13]) is presented in Section [H 
Section [5| is devoted to show our main result of estimating the rate of convergence (Theorem II. Sp . In Section [6l we 
give an application to the case of linear transport equations (Theorems ll.8l and ll.9p . We end up in Section [7| with 
an Appendix where we give the proof of Examples 11.31 and 11.61 

2 Ergodicity and preliminary facts 

In this section we present the proof of Propositions 11.11 and 11.41 We first start with the following ergodicity result 
which is a particular case of [121 Proposition 4.2]. However, we give the proof in our particular case for the sake of 
completeness. 

Proposition 2.1 (Ergodicity). Let g{v,T) : ^ M 6e a function satisfying: 
• (HI) Regularity: g is Lipschitz continuous with \\g\\i^^(jg2j < (3; 



4 



• (H2) Periodicity: g{v + 1,t + k) ^ g{v, r) for any (l, k) € , {v, r) G R^. 

Let V be the solution of the following equation 

jvr=g{v,T), r>0, 

1^(0) =.0, ^ ^ 

then there exist a constant A G M (independent of the initial data Vq) such that for every t, t' > 0, we have: 

\v{t) -v{t') ~ X{t -t')\ <^ with 1 + 2/3. (2.2) 

Remark 2.2 Under our assumptions, it is possible (see \12\ Theorem 1.5]) to show the existence of a hull function 
/i : M X R — > M satisfying: 

h{T + l,x) = h{T, x) 
h{T,x + 1) = h{T,x) + 1 
h, > 0, 

such that U{t, x) ~ h{T, x + It), with I given in Provosition \2. 1[ is a viscosity solution of: 

(Ur=g{U,T) 
\ U{t ^0,x)= x. 

This function h may he discontinuous, but its existence suggests that 

v{T) = h{T,lT) (2.3) 

is formally a classical of \2.1\) . The expression 112. 3\} allows to understand estimate 112. 2\) and also suggests that 
v{t) — It could be quasi-periodic in some cases. We emphasize the fact that in the proof of Provosition \2. 11 we do 
not use any notion of hull function, but propose a completely independent proof. 

Proof of Proposition 12.11 For any T > 0, define the two quantities: 

^+(r) = sup "(" + ^)-"(") and r(r)^inf <- + T)-vir) _ 

r>0 T ^ ' r>0 T 

These quantities are finite since Vt- = g is bounded. The proof is divided into three steps. 
Step 1: Estimate of \l+ - ^ |. 

Let 5 > be an arbitrary constant. From the definition of 1^{T), there always exists such that 



;±(T)- 



t;(r± +T) -t;(T±) 



< 6. (2.4) 



T 

Denote by [ij and [i] as the floor and the ceil integer parts of the real number t respectively. We consider 

fc=[r~— T+J, T~ = T~ — k and I = \v{t~) — v{t~)~\, (2-5) 
and consider w{t) = v(t + k) + I. Using (|2.ip and (H2), we check that is a solution of the following equation 

{Wr — g(w, t), r > 

^_ 1 (2.6) 

w{t ) = v{t ) + I. 

We remark, from (|2.5p . that 

T+ <T- <t+ + 1, (2.7) 
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and 

v{t-)<w(t~) <v{t^) + 1. (2.8) 
From (|2.ip . (|2.6p and (|2.8p . the comparison principle for ODEs gives in particular 

v{t)<w{t) for all t>t~. (2.9) 

If we suppose T > 1, we obtain from (|2.7p that t+ + T > and hence (using (|2.9p ): 

w(t+ +r) < u;(t+ +T). (2.10) 

Direct computations give 

{v{t- + T) - v{t-)) - {v{t+ + T) - w(t+)) = (u;(f" + T) - w(t+ + T)) + (w(f ") - w(?")) 

+ (u;(t+ + T) - w(r+ + T)) + (v{t+) - v{t-)), 

where, from ([SiHD and ((^117)) . we deduce that {v{t- + T) - v{t-)) - (u(t+ + T) - w(t+)) > -1 - 2||.g||oc- This 
inequality, together with ([2^ show that < /+(T) - /"(T) < and since this is true for any ^ > 0, 

we obtain 

l + 2||.g||oo 



\l+iT)-l-{T)\< 

However, in the case where T < 1, we always have \'"('^+'^^^'"('^'>\ < \\g\\r^ < Mhs. and therefore 

0</+(T)-r(T) < 



T 

T 

2||5l|oo 



T 

then 



|;+(T) - r(r)| < |; for every r>0. (2.11) 



Step 2: Existence of the limit of I'^iT) as T — > oo. 

First, if we compute 1+{PT) for P e N \ {0}, we get: 



Z+(Pr)-sup^ 

r>0 



>^ v{t + iT) - v{t + \)T) 

2-^ y 



Similarly, we get /~(PT) > /"(T). Consider Ti, Ta > such that TiP = T2Q for some P, Q G N \ {0}. Using this 
and (|2.1ip . we have 

1+{T2) > 1+{T2Q) = l+{TiP) > l-{TiP) > l-{Ti) > 1+{T,) - ^, 

J-i 

similarly we have Z+(T2) - ^+(Pi) < ^, then 



|Z+(Ti)-Z+(T2)|<max(l,l). (2.12) 



Pi P2 



By the same arguments as above we can get 



|r(ri)-r(r2)| <max(^l,^^ . (2.13) 

Recall that (|2.12p and (j2.13p are true when T1/T2 is rational. By an approximation argument, joint with the 

continuity of Z^, it is easy to see that this is still true when T1/T2 is any positive real number. Moreover, the 
identities (|2.12p and (|2.13p give that the sequence {1^{T))t is a Cauchy sequence as P — > 00, and hence it has a 
limit: 

lim 1^{T) = I, (2.14) 

T— >oo 
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which is the same hmit because of (P?TT|) . From (PH)) . (P?T^ and ^J^, inequaUty (|^ directly foUows. 
Step 3: Independence of vq. 

The fact that / is independent of vq foUows directly from the comparison principle and inequality (j2.2p . □ 
We can now present the proof of Propositions 11.11 and 11.41 



Proof of Proposition [TTTJ Using inequality (|2.2p of Proposition [5111 we can easily see that (|1.2p directly follows. 
It remains to show (|1.4p . We argue in two steps. 

Step 1: Monotonicity of 0/. 

Let u < u. Call I = 0/(u, t) and I — 0f{u, t). Let v and u be the solutions of: 

( Vr ^ f{v,T,u,t), T > 0, 
\i;(0) =iio, 

and 

\ Vr = fiv,T,U,t), T > 0, 
\ 5(0) = Uo, 

respectively. Assume without loss of generality that uq = 0. Using (A3), we deduce that f{v,T,u,t) < f{v,T,u^t). 
Hence, the comparison principle gives: 

v{t) < v{t) for every r > 0. (2-15) 
From inequality (|2.2p of Proposition [2Tll we have: 

\v{t) -It\<^ and \v{t) - It\ < ^ for aU r > 0. (2.16) 
We then easily conclude that / < Z as a consequence of (j2.15p and (|2.16p . 

Step 2: Continuity of 0/. 

We refer the reader to Proposition 13.21 which implies in particular the continuity of 0/. The main idea of the proof 

i. 

T- 



is to apply a perturbation argument using the inequality \l^(T) — l\ < ^- □ 



Proof of Proposition [TT4l 

Global existence. This is a direct consequence of the Cauchy-Peano theorem, using in particular the continuity of 
/ (see III- 

Uniqueness. Assume that there exists S C^([0, cxd); M) another solution of (|1.6p . Define k{t) = |u°(t) — u^(t)\, 
we compute (with the sign function sgn(a::) = a;/|a;| ii x ^ 0): 

hit) = {u",{t)~ulit))sgn{u''{t)^u\t)) 

= i0fiu'>it),t)~0f{u'it),t))sgniu'>{t)-u\t))<O, 

where for the last line we have used the monotonicity of 0/ (see (|1.4p ). This immediately implies that u° — u^. □ 

3 A stability result for the effective slope 0/ 

In this section, we will show a stability result for the term I given by Proposition 12.11 under a perturbation of (|2.ip 
of the form: 

j vt ^ gjiv,t) = g{v,t) ±j, t>0 
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where 7 > is a small real number. More precisely, we have the following proposition: 

Proposition 3.1 (Stability result). Take < 7 < 1. Let l-y be the effective slope given by Provosition \2.1[ which 
is associated to equation iS.l]) . and let Iq be the one corresponds to j = in iS.l]) . Then we have the following 
estimate: 

\lf-lo\<TT^ With e= (3 + 2e)(l + 2L), (3.2) 



and L 



Proof. We assume, for the sake of simplicity, that g^ = .9 + 7- In this case, it is easy to check that 1^ > Iq. The 
other case with = g — 7 is treated similarly. We first transform our ODE problem into a PDE one by setting v'' 
as the solution of the following equation: 

iv]{t,x)^g{v^{t,x),t)+-f, in (0,oo)xR 
\ v''{0,x) ^x, xeR. 

The proof is divided into three steps. 

Step 1: A control on wj. 

Using comparison principle arguments for (13. 3p . it is easily checked that v'^(t, .) is a non-decreasing function satis- 
fying v'^{t,x+ 1) = v''{t,x) + 1. We want to control v'J.{t, .) for any t. For this reason, we proceed as follows. Define 
for z > 0: 

r]{t,x) = v^{t,x + z) - ze^*, i > 0, a; G M. 

We compute 

r]t{t,x) = v]{t,x + z) - zLe^^ 

= g{r]{t,x) + ze^^t) - zLe^* + -f 
< 9{v{t,x),t) 

which proves that 77 is a sub-solution of (|3.3p with ri{0,x) ~ v'^(Q,x + z) — z — v''{0,x), and therefore, by the 
comparison principle, we obtain 

ri{t, x) = v'^{t, x + z)- ze^* < v'^it, x) 

hence for any t > 0, we have < v'*{t, x + z) — v'^it, x) < ze^*, then v'^ {t, x) is Lipschitz continuous in the variable 
X, satisfying: 

< v2it, x) < e^* for t > and a.c. x G M. (3.4) 
In a similar way, we can obtain a positive bound from below on vj, and finally get 

e-^' <v2ifx) <e^K (3.5) 

Step 2: An upper bound of v'^. 

We seek to find an upper bound of v'^ by constructing an explicit super-solution of p.3p with suitable initial data, 
and comparing it with v'^ . For this purpose, let 

w{t,x) ^v"{t,x + ci-ft), (t, x) e (0, 00) X K, (3.6) 

where vq is equal to v'^ for 7 = 0, and ci is positive constant to be precised later. We calculate: 

wt{t,x) = v'^{t,x + cijt) + Cl^V°{t,X + Cijt) 
= g{w{t,x),t) + cijv'^{t,x + ci'^t), 
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where from (|3.5p . we deduce that 

wt{t, x) > g{w{t, x),t)+ ci7e-^*. (3.7) 

Take ci = e^"^ for some fixed T > 0. Then using p.7p . we get wt{t, x) > g{w{t, x), + 7 for any t E [0, T]. Hence 
w is a super-solution of (|3.3p over [0,T] whose initial condition w{0,x) = v'^{0,x), which finally gives: 

w{t, x) > v^{t, x) Vt e [0, T], a; e R. (3.8) 

Step 3: Conclusion. 

We will now show the error estimate p. 21) . To this end, we will estimate both sides of inequality p.Sp involving Iq 
and l^. Firstly, using (|2.2p and (|3.5p . we compute: 

\v°{t, X + e^^7i) - v°(0, x) I < |w"(t, X + e^^7<) - v\t, x) \ + x) - v°(0, x) | 

We take this inequality for t — T and a; = 0, we get 

w{T, 0) = w°(T, e^^7T) < 7re2-^^ + /qT + ^. (3.9) 
Secondly, using similar arguments, and the fact that 7 < 1, wc obtain \v^{T, 0) — v'^{0, 0) — ljT\ < 2 + ^, hence 

v^{T,0)>l^T- {2 + 0- (3.10) 
Combining ^^E^, and (PrU]) . it follows that 

{l^ - la)T < jTe^^^ + 2{1 + ^) (3.11) 

Using p. lip , we deduce that: 

m-Zo|<7e^^^ + ^^. (3.12) 

Since the variable T was arbitrary chosen, let T satisfies jTe'^^'^ = 1 and therefore T > ^2£ll. From ([XT^ . the 
result directly follows. □ 

An immediate consequence of Proposition 13 . 1 1 is the following: 

Proposition 3.2 (Modulus of continuity of 0f). The function 0f{u,t) given by il.S]) satisfies for any {u,t) e 
and for all \v\ + \s\ < \: 

\0f(u^v,t + s)-0f{u,t)\<- ^ (3.13) 

|loga(|u| + |s|)| 

where a is given in assumption (Al). 

Remark that estimate p.l3p is optimal in view of Example 11.31 

4 Basic error estimate 

We start this section by considering a discrete scheme associated to the ODE (|1.6p . Namely, for a given w° (which 
may be chosen equal to uq or may be different), and for a time step At > 0, we define the sequence (w'^)fegN as 
follows: 

v''+^ = v'' + IkAt, lk = 0f{v^kAt), keN. (4.1) 

In this section we give a local error estimate between the solution u'^ of (jl.ip and the sequence v'^. To be more 
precise, we will show the following proposition: 
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Proposition 4.1 (Basic error estimate). Under assumptions (Al), (A2) and (A3) on the function f, let 
At > be small enough (depending only on a and (3), and take 

ei^\u%iM)-v% i = 0, 1. 

Then we have: 

cAt , , 

ei<eo + ce+- — , 4.2 

|logcAi| 

where c = c(q!, /3) > is a positive constant. 

The proof of the above proposition will be presented later in this section. In what follows in this section and in 
Section [51 we will assume that the arguments of the various logarithms are all less than 1 . 



Lemma 4.2 (Refined basic error estimate). Under the same hypothesis of Proposition pTTl let 

e+ = max(0,u'(iAt) - u*), e~ = min(0, u^(iAi) - d'), i = 0, 1 

and let = sup (max(0, u'^(t) — v^)), d(^ — inf (min(0, {/"^(t) — Then we have: 

te[o,At] te[o,At] 

e+ < e+ + (2 + Oe + . (4.3) 

and 

er > eo - (2 + e)e - > (4.4) 

\\og a{dj + At)\ 

with ^ = 1 + 2/? and [3 defined in (Al). 

Proof. In order to get estimates (|4.3I) . (14. 4p . the main idea is to freeze the last two arguments of /, and to use 
some comparison arguments. We start by estimating the term / |, u*^, t) from below. Since u'^{t) <v^ + d'l for 
t £ [0, Ai], we deduce, using (A3), that 

/(f,J,<t)-/(!f,i,.o,o)>/(f,J,.o + <,t)-/(f,^,.o,0 

and hence, from (Al), we get 

f(^^^u\t\>f(-iyA- a{dt + At). (4.5) 



Using similar arguments, we can also show 



e e 



We know from the definition of Cq and Cq that: 



f{-,-,u^t) <f(^,i,v0^0]+a{\d^\+At). (4.6) 



v° + <u'{0) <v" + e+. (4.7) 
Let and be the solutions of the following ODEs: 

= / ( — ,-,t;0,0 ) -a(d+ + At), t>0 



and 



e e y ^ ^ ' (4.8) 

vf{0)^v" + eQ 



wt^f(^^,-^,v°,0^+ai\d^\+At), t>0 
^v/{0) = v° + e+, 



(4.9) 
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respectively. From (j4.5p . (|4.6p and (|4.7p . we deduce, using the comparison principle, that: 

W'<u''<uf on [0,At]. (4.10) 

Applying Proposition 12. II to the functions W'^ and uf, we know that there exists two real numbers Iq and ^.q such 
that for t G [0, At] (assuming a(d^ + At) < 1 and a{\d^ \ + At) < 1): 

\vif{t)-uf{0)-lot\<{2 + Oe and \w'{t) - w'{0) - Igtl < {2 + ^)e. (4.11) 

Inequalities (|4.10p and (|4.1ip give: 

v° + e- +Jot-{2 + ^)e<u^{t)<v° + e+ +y+{2 + ^)e, te[0,At]. (4.12) 

Using Proposition 13. 11 we obtain, for Iq — 0f{v^,O), that 

0>Io-^o>-^ 7J^-X7^i (4.13) 

I loga(dj^ + At) I 

and 

0</n-^o< • (4.14) 

--° °- |loga(|dr| + At)| ^ ' 

The above two inequalities, together with (14. 7p and (I4.12p give the result. □ 

Two immediate corollaries of the above lemma are the following: 
Corollary 4.3 (Refined estimates involving ef). Under the same hypothesis of Lemma\J^.S\ we have: 



s+<e+ + (2 + 0e+ — (4.15) 

' - " ^ ^ |logai(|eo| + At)| ^ ^ 



and 

logai(e+ + At)| 

where ai = ai{a^[3) > is a positive constant. 



er>eo -(2 + Oe- '„+ , (4.16) 



Proof. We have d'^ < + (3At and \di | < Icq | + /3At (recall that (3 — ||/|1l°o(r<i)). Therefore, the result can be 
deduced from and ((i^ . □ 



Corollary 4.4 (Refined estimates with continuous time). Under the same hypothesis of Proposition /or 
every t G [0, At], define the continuous function e by: 

eit) =u'{t)~ {v° + lot). (4.17) 

Also define e'^it) = max(0,e(t)) and e^(t) — min(0,e(t)). Then we have for < ti < t2 < At and c = c{a,f3) > 0: 

e+(t2) < e+(ti) + ce+- ^f^ '^^l TT^ (4-18) 

|logc(|e (ti)| + (t2 -ti))| 

and 

e-(t2) > e-(ti) - C6 - , -T^. (4.19) 

|logc(e+(ti) + (t2 - ti))| 

Proof. We apply the same proof (word by word) of Lemma [4.21 and CoroUarv 14.31 with the origin shifted to ti 
and the time step At replaced by t2 — ti . □ 



Now we are ready to prove Proposition 14.1 
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Proof of Proposition 14.11 From the definition (|4.17p of e, we know that 

eo = |e(0)| and ei = \e{At)\. 
The case where ei = e(Ai) = is obvious. Four cases could be considered. 

Case 1. fe(0) > and e{At) > 0). In this case we have eo — Cq , ei — Ci and Cq = 0. Therefore (|4.2p is an 
immediate consequence of (|4.15p . 

Case 2. fe(0) < and e{At) <0). Similar to Case 1. 

Case 3. fe(0) < and e(At) > Oj. In this case ei = e+(At). Let the time t |_ be defined as follows: 

t_+ = max{i e (0, At); e{t) = 0}. 
Using inequality (|4.18p with ti = t ^, ^2 — At and At = At — t ^ < Ai, we get: 

cAt cAt 
ei < ce + T- — < Co + ce 



|logcAt| - |logcAt| 
where we have used the fact that e+(<i) = e~{ti) = 0, and hence (|4.2p follows. 

Case 4. fe(0) > and e(At) < Oj. Similar to Case 3. □ 

5 Estimate of the rate of convergence 

This section is entirely devoted to the proof of Theorem 11.51 Let T > and let At > be such that 

nAt = T, nen, (5.1) 

where n to be chosen large enough (the choice of At will be given later). In order to estimate \\u'^ — u^\\l°=(q,t)i '^^ 
add and subtract the continuous piecewise linear function passing through the points {kAt^v^), k = ■ ■ ■ n. In 
other words 

v{t) = v*' + {t - kAt)lk, kAt <t < {k + l)At, fc = 0---7i-l, 
where v'^ and Ik are defined in (|4.ip . We start by stating the following corollary that generalizes Proposition 14. II 

Corollary 5.1 (Basic error estimate in continuous time). Under assumptions (Al), (A2) and (A3), let 
T > 0, and At > given by i5.1\) . Define the function e{t) by: 

eit) ^ \u^t) ~ v{t)\, ie[0,T]. (5.2) 

For k — ■ ■ ■ n — I, call — e{kAt). Then for kAt < t < (k + l)At, we have: 

, , cAt 

et <efc + ce+- — , 5.3 

I log c At I 

where c — c(a,(3) > is a positive constant (the same given by Proposition \4-l\ l- 

Proof. We simply apply Proposition 14.11 with the origin shifted to kAt and with the time step At replaced by 
St^t- kAt < At. □ 

At this stage, we can show an inequality similar to (j5.3p with e(t) and replaced respectively by the functions 

eO(t) = \u°{t) ~ v{t)\ and el = e°{kAt), (5.4) 
where is the solution of (|1.6p . Indeed, we have the following proposition: 
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Proposition 5.2 (Basic error estimate for the homogenized ODE). Let 0/ be the function given by il.S]) . 
and enjoying the properties given by Provositions I j . il and \3.'A Let T > 0, and At > given by 115. Then for 
kAt <t<{k + l)At, fc = • • • n - 1, we have: 



e%t)<el + 



cAt 



logcAil 



where e", are given by |5.^[ j, and c = c(a,/?) > is a positive constant. 



(5.5) 



Sketch of the proof. Although the proof is an adaptation of the proof of inequahty (|5.3p (it suffices to deal with 
0f{u,t) instead of f{v,T,u,t), and to take e — 0), we will indicate the main points where it slightly differs. The 
crucial idea is that, on the one hand, the monotonicity of the function f(v,T,u,t) with respect to the variable u 
(see assumption (A3)) is replaced by the monotonicity of 0/(u,i) with respect to u (see (11.41) ). On the other hand, 
the fact that / is Lipschitz continuous (see assumption (Al)) is replaced by the fact that 0/ satisfies a modulus of 
continuity (see Proposition 13. 2p . and the fact that 0f{u,t) — l{u,t). 

Let us go into the details. In fact, the proof of Lemma [42] can be adapted where inequalities (|4.5[) and (|4.6p are 
replaced by 



and 



0/(w^^)>0/(^^^O)- 

0/(«°,t)<0/(«°,O) + 



— 0+ TTT' "^^^^ c^i'^ = sup (max(0,M°(t) - w^)) 

\oga(d{ + At)\ te[0At] 



with inf (mm(0.u°(t) -v°)) 

te[a,At] 



(5.6) 



(5.7) 



log«(|d;-| + At)r 

respectively. Here we have used the monotonicity of 0/, and inequality (j3.13p given by Proposition 13.21 Having 
(15. 6p and (|5.7p in hands, the sub- and super-solution w'^, uf defined by (|4.8p and (|4.9p are replaced by uP , 
solutions of 



and 



respectively, and we have 



with 



w° = 0/(«",o)- 

w°iO)^v° + e"^- 
w°=0f{v'>,O)+ 



loga(d?'+ + Ai)| 
with e° - = min(0, u"(0) - v°) 



loga{\4'-\+At)\ 
with e°'+ = max(0, u°(0) - v°) 



if{t)-w°{0)-llt = and w°{t)-w"{0)-l°t = 



and !^ — = lo + ■ 



loga(|d?'-|-t-Ai)r 



(5.8) 



(5.9) 



|loga(dy'+ + At)| 

where we recall the reader that Iq = 0/(w°, 0). Remark that (|5.8p replaces (|4.1ip . while (|5.9p gives (as a replacement 
of (jiJa and gH): 



> 7" - ^0 



loga(d?'+ + At)| 



and Q<Lo~lo 



\\oga{\d',-\ + At)\- 



At this point, the rest of the proof proceeds in a very similar way as the proof of (j5.3p with e = 0, and a possible 
changing of the constants but always depending on a and /?. □ 



Now we are ready to present the proof of our main result. 
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Proof of Theorem 11.51 We first note that the constant c > may certainly differ from hne to line in the proof. 
We decompose the quantity \u'^{t) — u"{t)\ into two pieces: 

e(t) e«{t) 

\u'{t) - < '\u'{t)-v{t)\ + \u"{t) - V{t)\ . (5.10) 

In order to estimate e(t), we iterate inequality (|5.3p and we finally obtain for m'^(O) = uq = v'^ , and t e [0, T] with 
T = nAt: 

cnAt ceT cT 

e(t) < en + cne + t- r— r < — ; h -n 1— r- 

^ ^ - |logcAi| - At |logcAt| 

The above inequality gives: 

and, from inequality (|5.5p of Proposition 15.21 we can show in the same way as above that we also have: 

Choosing particularly At = Ce| loge|, we deduce that e < At = Ce| log e| < T, and (from (|?TTT|) . ([5?T^ ) that: 

e(t)<^^ and e°(t) < (5.13) 
|loge| |loge| 

with c in (j5.13|) depending on the choice of C > 0. finally, inequality ()1.7|) could now be easily deduced from (|5.10p 
and ([5l3l) . □ 

6 Application: error estimate for linear transport equations 

In this section, as an application of our previous results on ODEs, we give the proof of some error estimates for the 
homogenization of linear transport equations. Namely, we prove Theorems 11.81 and II .91 We start with Theorem ll.8l 
keeping the same notations of Subsection 11.21 



Proof of Theorem 11.81 The proof is divided into four steps. The first three steps are devoted to the definition 
of the limit solution V'^, and to prove that it is a viscosity solution. The proof of the error estimate is done in the 
last step. 

Step 1: Definition of V°. 

Because the homogenized vector field 0a is not Lipschitz, we define our solution to (II. lip in an indirect way 
using the characteristics. Precisely, for {t,x) G (0,oo) x K.^, x = {xi,X2), we define V'^{t,x) as follows: 

y°(t,x)-yo(XO(0;t,a:)) (6.1) 

where the curve X°{t; t, x) : t e M ^ X°{t; t, x) e R^, 

X°iT;t,x) = (XO(T;t,x),XO(r;t,a:)) (6.2) 

is the solution of the following ODE: 

r a.X" = 0a(X") 

\xOit)^x. ^'-'^ 

We will see below that this solution is unique. For the sake of simplicity of notations, we will omit the dependence 
of X'^ on (t, x) and we will simply write 

X^{T;t,x)^X"{T). 
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From (|6.3p and (|1.9p . we can easily check that ^2 (''') — ^2 ^ t + t, hence using (|6.ip . we get 

= t/o(^i"(0),X2-i) 

where Xj' satisfies: 

r a,X? = -0/ (X°,a;2-t + T) 

In order to show that X^ (0) is uniquely defined, we solve (|6.4p backwards, in other words, we let 

0X°(r) = X°(t-r). 



(6.4) 



In this case: 

where 0X'^ satisfies: 



V''{t,x) = Vo(0X°,{t),X2-t), (6.5) 



dr0X° = 0f {0X°,X2-t) 

0X°{O)^xi. 



(6.6) 



From Proposition 11.41 the solution 0X'^ e C^([0,oo);R) is unique and hence X^{0) = 0X^{t) is uniquely deter- 
mined. Consequently the function V'^ is well defined. 

Step 2: is Lipschitz continuous. 

From Step 1, we know that 

V''{t,xi,X2) ^Vo{0X^{t),X2-t) (6.7) 

with 0X1 given by (|6.6p also depends on xi and X2- Let the function y : M'^ ^ M be defined as follows (with 
simplified notation showing the dependence on the variables (t, xi, X2))- 

Y{t,xi,X2) :=0XO(i). 

In order to show that y° is Lipschitz, it suffices (see (|6.7p ) to show that Y is Lipschitz. First, it is easily seen from 
(|6.6p that Y is Lipschitz in time t. The Lipschitz continuity with respect to the variable Xi directly follows from 
the monotonicity of 0/ (see (|1.4p ). and the comparison principle. In order to show the Lipschitz continuity with 
respect to X2, we first give a formal proof by assuming that 0/ is smooth, and then we present the main idea that 
permit to make the proof rigorous. Suppose that 

0/eC°°(M2;M) with 10/(2/, s)| < II0/IIOO. 

Take 



then the above two functions satisfy: 



and 



Y = da^^Y and Y = dtY, 

(dtY^{dy0f)Y + d,0f ^^^^ 

\y{0,xi,X2) = 0, 



dtY = idy0f)Y ~ ds0f 

(6.9) 

F(0,a::i,a;2) = 0f{xi,X2), 



respectively. Let 0Y = Y + Y, we get (from (|6.8p and (|6.9p ): 



f dt0Y = {dy0f)0Y 

I 0y(O,a;i,a;2) = 0/(a;i,a;2), 
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which gives, because of the monotonicity of 0/ (see (|1.4p ). that the function t \0Y{t, ., .)| is non-increasing. Hence 

\0Yit,xuX2)\ < \0fixi,X2)\ < II0/IU. (6.10) 
Since Y = dtY, we know from (|6.6p that |y| < |l0/||cx; where we finaUy obtain (see (IG.lOp ): 

\Y\^\d^,Y\<2\\0f\\^, 

which shows that Y is Lipschitz continuous in the X2 variable. In order to make the proof rigorous, it suffices 
to consider a regular approximation of the function 0/ (as for example the convolution with a suitable moUifier 
sequence) and then to pass to the limit. 

Step 3: V'^ is a viscosity solution of (jl.lip . 

For the definition and the study of the theory of viscosity solutions, we refer the reader to |3]. Let us take 
0, e C^(IR^;K) such that V''^ — 0(j) (resp. V'^ — (p) has a local maximum (resp. local minimum) at some point 
(0t,0x) G (0, cxd) X 'M? (resp. {t,x) G (0, 00) x R^), with V°{0t,0x) = 04){0t,0x) and V°{t,x) = (jAt^x). In order to 
show that is a viscosity solution of (jl.lip . we need to show the following two inequalities: 

dt04'{0t, 0x) + {0a ■ V0(/))(0t, 0x) < (6-11) 

and 

dt±{hx) + (0a • V0)(i,x) > 0. (6.12) 

We only show inequality (|6.1ip . In fact, inequality (|6.12p can be proved in exactly the same way. For any t G [0, 0t] 
define the function (j) by 

</> : t -> (l){t) = 00(t, X"(t; 0t, 0x)), 

where X° is defined in (|6.3p . Let us show an inequality on in the interval [0t — r,0t] for r > small enough. 
Remark that X'~'{0t; 0t, 0x) = 0x. Hence, for t E [0t — r, 0t], (t, X^{t; 0t, 0x)) is close to {0t, 0x) and therefore (since 
V'^ — 0(j) has a local maximum at {0t, 0x) with T^°(0t, 0x) — 0(j){0t, 0x)) we get: 

0(i) = 00(t, X°{t; 0t, 0x)) > y"(t, X°{t; 0t, 0x)) = Va{X°{0; t, ^^(i; 0t, 0x))) 
= Vo{X"{0; 0t, 0x)) = V^i0t, 0x) = 0(t){0t, 0x) = 0(t){0t, X^{0t; 0t, 0x)) = 4>{0t), 

where the passage from the first to the second line is due to the fact that the points {t,X^{t\0t,0x)) and (0t, 0a:) 
are on the same characteristics. Finally, this implies 

which directly gives (|6.1ip . 

Step 4: Proof of the error estimate (|1.13p . 

The solution of (jl.lOp can be written (in analogue with ()6.5p and ()6.6p ) as 

V'{t,x) = Va{0Xl{t),X2~t) (6.13) 

where the characteristics 0X{ satisfies: 

a,.x:.;(?S,-_I,X..,.,_.) ^^^^^ 

0Xl{Q)^xi. 

We apply Theorem 11.51 namely inequality (|1.7p . with u*^ and replaced by 0XI and 0X^ respectivly, we obtain: 

||0Xi^-0X°|L < for r>Ce|loge| with C>0,e>0. (6.15) 

" ^^'-^1 |loge| 
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Using (lEni), (Pl^ and (HHD, we compute for {t,x) e (0, T) 



p2. 



\V'{t,x)-V"{t,x)\ = |T/o(0^i'W,a;2-t)-'K)(0X?(t),a;2-t)| 
< Ltp{Vo)\0Xl{t) ^ 0X^(t)\ 

|loge| 

and inequality (|1.13p directly follows. 

Proof of Theorem II. 9i Application of Proposition [27T1 and same proof as Theorem [T 

7 Appendix: proof of Examples 11.31 and 11.61 

Proof of Example 11.31 The function 0/ can be expressed as (see for instance |23|): 



0f{u,t) 



dv 



Q — u + I sin27rw| 



f ^ dv f dv 

Let a = —u > 0, it is easy to check that / j — 2 j -. Take 

Jq a+|sm27ru| 7|„|<i/4 a + | sm 27rw| 

r = / . and /-^ ^ 



l\v\<i/4 a + I sin27ru| J\v\<Ra a + 2n\v\ 

We are interested in the limit a — > and R ^ oo with Ra — > 0. We compute 

A B 



ja _ ja,R ^ f dv ^ r 27r|t>| - |sin27rt;| ^^^^ 

yi?,a<|t,|<i/4 a+ |sin27rw| y|„|<fl„ (a + 27r|'(;|)(a + | sin27rw|) 

where we have: 

B ^0 as Ra^O 
11 c 



A < — c\/\ logal 

- 2 sin(27ri?a) Ra v i b i 



with c = i^- and we have chosen 

R ' 



i\/|loga| 
Now, let w = -, we also compute: 



/"•^ = 2 / ^^ = i(log(l + 27ri?))~llogi?^l|loga| 
Jo<v<R 1 + Sttw tt tt tt 



/0<v<R 

Since A < /°'^, this shows that I" ~ /°'^ and then 

^ oil ^ I II as M ^ 0", 
2|log|w|| 

which justifies (|1.5p . 

Proof of Example 11.61 Since /(u, r, u, t) ~ g{v + r) — 1, then the function u*^ defined by: 
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satisfies 



vl^g[-], t>0, 



ej' (7.1) 

v'{0) = 0. 

From the particular expression (|1.8[) of g with g (^) — 0, we can check that < ^ < i for every t > 0, which 
imphes that 

v'-^O in 



then 
where 

Moreover, equation (|7.ip can be written: 



u° in L°° with = -t, 



1 

= with u^(0)=0, 

2 e w > 



therefore (solving the above equation) we get for t = Se] loge| 



^it)-u\t)=v^it) = |(l-e-^) 

2^ 



2 2 2 25|loge|' 

which terminates the proof of Example 11.61 □ 
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